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TECHNICAL MEMORANDUM X-64541 


CALCULATION OF THE PROFILES OF THE STOKES 
PARAMETERS AND THE DEGREE OF LINE POLARIZATION: 
AN APPLICATION OF THE SOLUTIONS OF MOE TO THE 
UNNO TRANSFER EQUATIONS 


SUMMARY 


This report outlines a particular solution to the Unno transfer 
equations and describes in detail its application to the calculations of the 
line profiles of the Stokes parameters and to the calculations of the linear 
and circular polarization integrated over the bandpass of a birefringent 
filter. 


Unno has derived the four transfer equations for the Stokes 
parameters that describe the state of polarization of a solar absorption line 
formed in a homogeneous magnetic field. Moe has found a solution to these 
equations under the assumption that the normalized function describing the 
variation with wavelength of the line absorption coefficient is independent of 
optical depth. These solutions must be numerically evaluated using a model 
of the solar atmosphere that gives the variations of temperature, gas 
pressure, and electron pressure with optical depth. 

In the numerical calculation, three depth-dependent functions are 
calculated: the source functions, the continuous absorption coefficient at 

Q 

the wavelength of the absorption line , and the line absorption coeffi- 
cient at line center with zero damping (0) . The line is assumed to be 

formed by the mechanism of pure absorption so that the source functions 
become identical to the Planck distribution function which can be written 

with an explicit temperature dependence. K is calculated in a sub- 

A 

program whose form is dictated by the model atmosphere. The "fictitious" 
line absorption coefficient K* (0) is calculated by assuming that the 

A 



population density of the lower level of the absorption transition can be 
calculated from the Saha and Boltzmann equations; this implies the assump- 
tion of local thermodynamic equilibrium. 

The wavelength-dependent part of the line absorption coefficient, 
which Moe has assumed to be independent of optical depth, is calculated by 
assuming contributions from natural, collisional, and thermal broadening 
mechanisms. It is assumed that these can be treated independently so that 
the resulting profile is the Voigt function which has been numerically 
tabulated for different values of the total damping parameter. 

The linear and circular polarization at a particular wavelength in 
the line is obtained from the linear (Q) and circular (V) Stokes parameters. 
In a magnetograph system using a narrow-band filter, the polarization that 
is measured is the contribution from all wavelengths in the interval of the 
bandpass. The actual measured quantities are thus numerically evaluated 
from the solutions for the Stokes parameters (as described above) and from 
the transmission characteristics of the filter. 

Numerical results are calculated for the Stokes parameters ?ind the 
degree of linear and circular polarized light for the case of the iron line 
5250. 22 A using the photospheric model of Holweger. 


INTRODUCTION 


Recently Moe [1] has obtained a solution to the Unno [2] transfer 
equations for the Stokes parameters, assuming only that the magnetic field 
is homogeneous and that the line absorption coefficient profile is independent 
of depth. In this paper the applications of this solution are described first 
for the calculations of the line profiles of the Stokes parameters and, 
secondly, for the calculations of the degrees of linear and circular polariza- 
tions that would be theoretically measured by a real-time solar magneto- 
graph system using a narrow-band birefringent filter. 


EQUATIONS DERIVED BY MOE 


The transfer equations of Unno [2] as modified by Moe [1] to 
include contributions to the Stokes parameters from scattering as well as 
from thermal emission are: 
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(1) 


In these equations I, Q, and V are the Stokes parameters that describe 

the state of polarization of light of intensity I (in the Unno coordinate 

system, for a homogeneous magnetic field, the fourth parameter U is 

zero) ; ju = cos 6 where 0 is the angle between the solar radius and the 

direction of the light beam; t is the optical depth in the continuum; and 

c 

S , , and S are the three source functions for the Stokes parameters, 

i y v 

The functions rj , rj , and rj are defined as follows: 
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sin 2 ip 
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COS ip 


and 


(2) 


In these equations, ip is the angle between the magnetic field vector and the 
c 

line of sight, K is the continuous absorption coefficient at the wavelength 
A 

of the line, and , and are the absorption coefficients for 

linearly, right-circularly, and left-circularly polarized light respectively. 


3 



It is assumed that these coefficients are identical to the line absorption 

coefficient k| (A) in the absence of a magnetic field except for the Zeeman 

displacement AA in wavelength; for a simple Zeeman triplet, 
xi 


± AA = 4. 67 x 1 0~ 5 gA 2 H 
H ° 


(in centimeters with A 0 in centimeters) , where g is the Lande g factor 
and H is the strength (in gauss) of the magnetic field. Thus, 

KMA) = Kj(A ± AA^) = K f (A+ AA fl ) = K^(A) 

The line absorption coefficient in the absence of a magnetic field, K* (A) , 
can be written as 


k} (A) - kJ (0) H(a , v) 

A A 

where (0) is the line absorption coefficient in the line center for zero 

damping, and H(a , v) is the normalized shape function for the line 
absorption coefficient. H(a , v) depends on the damping parameter "a" 
and the variable v= (A - A 0 )/AA^ , where AA^ is the Doppler (or thermal) 

broadening factor, 


» -v An / 2kT 9 
D c n M s 


(3) 


In this expression, c is the speed of light, k is the Boltzmann constant, 

T is the absolute temperature, M is the mass of the atom being considered, 
and 4 is the microturbulent velocity of individual elements of the emitting 
gas. The function H(a , v) is usually assumed to be a folding of a 
Lorentzian (damping) profile and a Doppler profile, 


+ oo 


-v 2 

H / a - a r ej_dy__ 

H(a ’ v) ~ 7T J a 2 + (v - y) 2 
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various authors have described this function. (See Jefferies [3] for a list 
of references. ) 


From the above discussion, the following can now be written: 


K k:(o) 

7T X 

c c 
K x *X 


H(a, v) 


K r K X <0) 


H (a , v ± 


K l V°> 


H a , v + 


Substituting equations (4) into equations (2) and defining V = — — , 

H AA D 

the following are obtained for the Unno variables tj , t? and n : 

i y v 
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In his solution, Moe assumes that 


< ( °> 


K. 


is independent of wavelength 


(which is valid since the continuous absorption varies extremely slowly 
over the wavelength interval of the absorption line) but dependent on optical 
depth r , whereas F, 


I,Q,V 

dependent. The assumption that F 


are independent of depth but are wavelength- 
do not vary with depth imposes 


i.Q.v 

the constraint that the magnetic field must not vary with depth; i. e. , it must 
be homogeneous. It also means that H(a, v) must also be constant with 


depth, a condition that implies the damping ("a") and thermal 


(AX d in v) 


broadening parameters cannot vary with depth in the expression for H(a,v) . 
Under these assumptions, Moe finds the following solutions for the emergent 
intensities of the Stokes parameters: 


I =-(Z 1+ Z 2 ) 


and 
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APPLICATION TO THE CASE OF 
PURE ABSORPTION 

If it is assumed that the Fraunhofer line is formed by the mechanism 
of pure absorption. 


S I - S Q - S V ' V T) 

where B (T) is the Planck intensity distribution function. Under this 
A 

assumption, equations (6) and (7) reduce to the following: 
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c 1 

Inequation (8), the depth-dependent functions are K , K (0) , and 

A A 

B(t ) ; the wavelength (and magnetic field) dependency is found in the 
c 

functions Fj y . To compute Z t and Z 2 , and thus the Stokes 

c 1 

parameters, the values of K. , K_ (0) , and B(r ) must be calculated at 

A A C 

each depth in the solar atmosphere; this means a "model atmosphere" 
must be employed that gives the temperature (T or 6 = 5040/T) , gas pres- 
sure (P in dynes per square centimeter) , electron pressure (P in dynes 
g e 

per square centimeter) , and continuous absorption coefficient at a standard 
c 

wavelength (K Q ) at various increments of optical depth . Since r ^ is 

a function of wavelength, model atmospheres are calculated for a standard 
wavelength; to change the depth scale to the wavelength being considered, the 
relation 


dr 

c 

O I 
T3 


w 

o o 


is used, where represents the continuum optical depth at the standard 
wavelength. The last of equations (8) then becomes 
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Finally, for computational purposes, the variable 


(9) 


= log T r 


is used rather than . In this notation the integrals in equation ( 9) 
transform into 



jn Mod 


where Mod=log 1 ^e . 

The intensities of the Stokes parameters are usually expressed 

relative to the background continuum. The expression for the continuum 

intensity I can be obtained by letting F _ =0 inequations (10) 

c 1, y, V 

and setting I = ^ (Z + Z ) ; the result is 

C A 1 Z 
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Equations (10) and (11) are the basic formulas from which the 
line profiles of the Stokes parameters are calculated. Of the terms in 
cl 

these equations, K , K (0), and B(x ) must be calculated from their 
A A vl 

Q 

known dependence on 6 , P g , and P g ; K Q (x 0 ) is (usually) given with 
the model. 


CALCULATION OF THE PLANCK FUNCTION 

The Planck function is given by 


V T > - 


2hc 


hc/AkT , 
e - 1 


Because our final expressions will be normalized by the continuum intensity 

2 


I , the constant factor 
c 


2hc 


can be ignored and the following written for 


B, 


( T( V) : 


B( V 


1 

(2.85476 x 10 4 )( s a 0 ) 
e - 1 


> 


with (the central wavelength of the absorption line) expressed in 
angstroms. 
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CALCULATION OF K° x 


In the construction of a model atmosphere, the continuous absorption 

Q 

coefficient at the standard wavelength, K Q , is either directly calculated by 

assuming a certain number of sources of opacity (such as the negative 
hydrogen ion H~ , neutral hydrogen H , metals, etc. ) or is obtained 

Q 

from tables of K Q already calculated by other authors (e.g. , Vitense [4], 

Bode [5] ) . In the first case, a subprogram must be constructed that 

c c 

calculates K using the same opacity sources used in calculating K 

A 0 

Q 

(and that calculates values of in agreement with those of the model) . 

In the second case, the same tables must be used that were used by the 
author of the model for the wavelength of the absorption line under study. 

In the actual calculations described in this paper, one of two methods can 
be used. The first method calculates from basic physical principles the 
continuous absorption coefficient with H and H~ as the sources of opacity; 
it has been found that by varying the helium to hydrogen ratio (B) in this 

Q 

subprogram, the computed values of can be made to agree with those 

of the model in the optical depth regions of interest. This subprogram has 
been described in detail by the author elsewhere [6]. The second method 
is a double interpolation (0 , P ) subprogram for the data of Bode [5] at 

G 

X 5250 and is described elsewhere by Wadsworth [7] . 


CALCULATION OF K.(0) 

A, 


x 


The line absorption coefficient K^(A) per gram of stellar matter, 
corrected for negative absorptions , is given by 


- N f B lu bv 0 


1 - e 


- hc/A Q kT \ 


/0(A) ; 


( 12 ) 
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here 0 (X) is the normalized profile of the line absorption coefficient. If 
a combination of natural, collisional, and thermal broadening is assumed, 


0(A) 





H (a , v) 


(13) 


For line center, zero-damping 


H ( 0 , 0) = 1 


and equations (12) and (13) combine to give 


. / -hc/\ kT\ X n 2 

< ( 0 »= N f B lu to «l 1 - e /-p 2 — 

NTT cAX 


D 


(14) 


or 

In equation (14) , is the number of atoms in the lower (1) level of the 

the line transition per gram (g) of stellar matter, is the Einstein 

coefficient of absorption, and hv . represents the energy difference between 

the lower (1) and upper (u) levels. The Einstein coefficient is related to 
the oscillator strength f (sometimes called the f -value) by the equation 


% 
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Using the definition for , equation (3) , and the notation 
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equation (14) becomes 
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(STEM) 


(15) 


K >> 


J 7 e \ f 


m c 
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+ I 
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g 

The number of atoms per gram of stellar matter, , can be expressed in 
terms of the number of atoms in the lower level per cubic centimeter, , 

N , 

N® = 

1 

where m is the atomic mass of hydrogen, e is the mean molecular 
H H 

weight of stellar matter relative to hydrogen, and N is the number of 

ri 

hydrogen nuclei (atoms, ions) per cubic centimeter. Regrouping the terms 
in equation (15) gives 


K ‘(°) 





(STEM) 


(16) 


To establish the population density of the lower level, local 

thermodynamic equilibrium is assumed and the Boltzmann and Saha equations 
are used. If it is assumed that 


N , = total number of nuclei of the element per cubic centimeter, 

El 

= the number of atoms of the element per cubic centimeter, 

N. = the number of singly- ionized atoms of the element per 
cubic centimeter, 

N.. = the number of doubly- ionized atoms of the element per 

cubic centimeter, etc. , 
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the following can be written (neglecting molecular structure): 



Also, the ratio N ,/N can be written as 

I H 


V n h ■ 



5 

N. 


N. / N.. N... N.. \ 

. 1 / . xi in n \ 

1 + n“ (‘ + n7 + 57 5T + - ' 7 

A v i li l ' 


N. 

N 


El 
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( 17 ) 


The first term in brackets in equation (17) can be calculated from the 
Boltzmann relation. 


Nj gj ~ X l/kT 



g^ is the statistical weight of the lower level, is the partition function 

of the atom, and x, is the excitation potential of the level. If N /N is 
i 1 A 

designated as N , then it can be calculated from the formula 


14 



-eXj/Mod 

Ni = — e 

A 


where x^ is in electron-volts. 


The second term in brackets in equation (17) is calculated by 
repeated applications of the Saha equation, 


N 


q + 1 
N 


N = 
e 


2 tt m kT\ 3//2 2 U 

e \ 9Ji± 


X q /kT 


U 


q \ h / q 

N is the number of atoms in the qth stage of ionization, N 


q + 1 


is the 


number in the (q + 1) stage, x is the ionization energy from q to 

(q + 1) , N is the electron number density, and the U's are the respective 
6 

partition functions. If it is assumed that N.. ,N... ,...<<N. (a valid 

li in l 

assumption for sunspot and photosphere models for most elements) , this 

second term reduces to 1/N where 

£ 


N. 

N 0 = 1 + -7- 

2 N a 


2U i T 5/2 -(01* + 0.48) /Mod 

1 + IT — 6 

A e 


Any effects of reduction of the ionization potential have been included here, 


I* = I - 4. 98 x 10~ 4 0 \Tp 


with P equal to the electron pressure in dynes per square centimeter and 

G 

"I" equal to the first ionization potential of the element in electron-volts. 

The third term in equation (17) is given from abundance determi- 
nations ( a e j) ; the assumed abundances used in constructing the model 

atmosphere should be used here. 
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Finally, equation (16) is written in the form 




where 





(STEM) 


K 


n/T 


m m TT c 
e H 



(18) 


(19) 


The "constant" K can be calculated directly for a given line or it can be 
determined by fitting the calculated profile for zero magnetic field to an 
experimental profile; the latter procedure has been adopted. 

In the calculations of N and N , the partition functions U . 

1 Z A 

and U. must be calculated as a function of temperature and electron 

pressure. This is done in this paper by a subprogram that evaluates the 
partition functions following the method of Traving et al. [8]. This sub- 
program was adapted from a program given to the author by Olav K. Moe 
(private communication). 


THE MODEL ATMOSPHERE 

To evaluate the integrals in equations (10) and (11) , it is 

cl 

necessary to be able to calculate , K^(0) , and B(x^) at discrete 

Q 

intervals of x„ ; this means that values of 0 and P (and KJ must 
0 e 0 

be obtained at these intervals. To use any desired increment in x Q , least 

c 

squares polynomial fits to the 0 (x ) , P (x ) , and K (x ) data given 

U 6U U U 

in various model atmospheres have been made, and these results have been 
incorporated in three subroutines into the main computer program. 
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THE FUNCTIONS F 


l,Q,V 


The remaining terms in equations (10) and (11) that must be 
evaluated numerically are the wavelength-dependent functions F _ 

defined by equations (5). These functions are determined from the magnetic 

field inclination ip and from the functions H(a,v) , H(a,v+v ) , and 

H 

H(a, v - v ) , where v = (A - A )/AA and v = AA /AA with A A 
ri U JJ H ri JJ D 

having a constant value. The Zeeman splitting factor AA is calculated 

from the assumed magnetic field intensity, and the Voigt functions H (a , v) 
are calculated from least squares polynomial fits to data tabulated by 
Hummer [9] for various values of the damping parameter "a". 


CALCULATED STOKES PROFILES 


Typical results are shown in Figures 1, 2, 3, and 4 for the Stokes 
parameters for the iron line 5250. 22 A . The constant K was given the 
value log K = 7. 914 , and the damping parameter "a" was 0. 10 ; the 
photospheric model of Holweger [10] was used. In this particular calcula- 
tion a variable velocity field £ was used in equation (18); Figure 5 shows 
this velocity field together with the velocity distribution calculated in 
Holweger's model. A constant Doppler half-width of 34. 0 mA was used in 
the wavelength-dependent functions. The measured profile in Figure 1 was 
obtained by Dr. John Harvey and the author using the Kitt Peak solar 
telescope and spectrograph in the photoelectric double-pass mode. 

In Figure 2, it is interesting to estimate the inclination angle ip 
using the Seares [11] formula, 

r : r =4(1+ cos 2 ip) : 4 sin 2 ip , 

a ir 4 r 2 r 


where r = 1 - I/I . From Figure 2, r = 1 - 0.6626 and r = 1 - 0.64 ; 
c a 7T 


0.338 


7 r 


0,36 


l • 2 

— am 6 


sin 4 ip 


This gives ip = 56. 5 deg compared to the actual value of 60 deg. 
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Figure i. Iron 5250. 22 A line p: 







Figure 4. Iron 5250. 22 





CALCULATION OF THE DEGREE OF LINEAR 
AND CIRCULAR POLARIZATION USING 
A NARROW-BAND BIREFRINGENT FILTER 


A solar magnetograph system, by the successive (or simultaneous) 
use of various polarizing optics and with a subsequent computer reduction of 

Q/I 


the resulting data, can provide measurements of the quantities 


I/I 


and 


i/i 


, averaged over the bandpass of the light-filtering system used. 


The transmission characteristics of the filter can be expressed by a function 
T (AA , AA. ) , where AA. represents A. - A Q with A. being the wavelength 

at which maximum filter transmission occurs. If the total bandpass of the 
filtering system is 2A , the measured quantities can be mathematically 
expressed as follows: 


AA. + A 

7 a( - f I & zf AX> ' T(*X . *,) d<*X, 

AA. -Ac v ' 

p = i 

Q, V AA. +A 

.. 7 . Tw T ( a * ■ ) d(ax) 

AA. -Ac ' ' 

l 


( 20 ) 

The integration is generally performed over the variable v = AA/AA Q ; 
defining 6 = A/AA^ and v. = AA./AA^ , equations (20) become 



v. + 6 

i 



v. + 6 

l 
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v. - 6 
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Q(v) t V (v) 
I c (v) 



I(v) 

! c ( v) 



( 21 ) 


23 



The interest here is in calculating 



for a solar magnetograph 


system using a Zeiss narrow-band birefringent filter. The theoretical 
transmission profile of such a filter is given by the equation 


(L-l) 

T(AX') = T ( 0) n 


k=0 



(22) 


In equation (22) , T(0) is the peak transmission, L is the number of filter 
elements, d/2 is that value of AX' where the profile has its first minimum, 
and Ax' is the wavelength distance measured from the peak transmission 
wavelength. Since all wavelengths (AX) are measured relative to the center 
(X Q ) of the absorption line, AX,' is written in terms of AX as follows: . 


Ax' = AX - AX. 

l 


If the variable v is introduced, equation (22) becomes 


T 




. (L-l) 

v. ) = T(0) n 
l ' k=0 



(23) 


The actual transmission profile of the Zeiss filter that is to be used 
in the NRL/MSFC real-time solar magnetograph was measured in December, 
1968 by J. Zwiener (private communication) using the Kitt Peak solar 
telescope and spectrograph; the data were reduced by W. C. Snoddy (private 
communication) . The profile at the central portion of the field of view was 
found to be slightly asymmetric with a bandpass of 0. 122 A. The asymmetry 
can be simulated in the theoretical profile by introducing a phase shift into 
the first term in the product as follows: 


T 



= T(0) 



l 

Q 

<1 

/ 

N 

fi 

COS 7T 

[(v- V ij d + «| 
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x 


(L-l) 

n 

k=l 


cos 


7T 





(24) 
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It was found that a good fit to the measured profile could be obtained for 
L = 8 , a = 0. 85 , and d = 0. 2752 A ; Figure 6 shows this comparison. 
The introduction of a phase shift a into the formula will shift the peak 
transmission wavelength from the value v by an amount depending on the 

magnitude of a ; the true filter maximum position is obtained by adding a 
correction term to v. . For a = 0. 85 , the shift is 0. 030 A and the 

l 

correction term is 30. 0/AA.^ . 

NUMERICAL RESULTS 


Calculations have been made using equations (21) and (24) for the 
degree of linear and circular polarization using the same line and parameters 
as in Figures 1 through 5. One particular application of these equations is 
in "tuning" the filter across the line; the filter's maximum position v is 


varied through the wavelength region of interest at a few values of H and ip. 

These results are of interest in determining optimum filter settings for 

linear and circular polarization measurements. A typical result is plotted 

in Figure 7. A simple modification allows the calculation of P at one 

Q 5 V 

or two particular filter settings v. for a large number of values of H and 

ip. These results are useful in "calibrating" magnetographs; calculated 
values of P are plotted versus H at different inclination angles ip to 

relate the magnetic field strength and direction on the sun to the measured 


quantities 


P and 

y 


V 


Figure 8 indicates a typical calibration curve. 
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PERCENT TRANSMISSION, T/T 







H (gauss) 


Figure 8. Magnetograph calibration curve, birefringent filter 
(filter setting at AX = 0. 051 A, 

inclination angle = 0 deg) . 
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